Gain, directionality and noise in microwave SQUID amplifiers: Input-output approach 
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We present a new theoretical framework to analyze microwave amplifiers based on the dc SQUID. 
Our analysis applies input-output theory generalized for Josephson junction devices biased in the 
running state. Using this approach we express the high frequency dynamics of the SQUID as a 
scattering between the participating modes. This enables us to elucidate the inherently nonreciprocal 
nature of gain as a function of bias current and input frequency. This method can, in principle, 
accommodate an arbitrary number of Josephson harmonics generated in the running state of the 
junction. We report detailed calculations taking into account the first few harmonics that provide 
simple semi- quantitative results showing a degradation of gain, directionality and noise of the device 
as a function of increasing signal frequency. We also discuss the fundamental limits on device 
performance and applications of this formalism to real devices. 

PACS numbers: 85.25.Cp, 74.50.+r, 85.25.Dq, 85.25.-j 



I. INTRODUCTION 

For almost half a century, the dc SQUID (Supercon- 
ducting QUantum Interference Device) has enabled a 
broad range of devices, including magnetometers, gra- 
diometers, voltmeters, susceptometers and amplifiers 
[1, 2]. Most of these devices are used at relatively low 
frequencies, and all have the common feature of offering 
extremely low noise [3]. The fact that the dc SQUID 
is potentially a quantum-limited amplifier in the mi- 
crowave regime was recognized long ago [4], but not ex- 
ploited in practice until the Axion Dark Matter experi- 
ment (ADMX) provided a powerful motivation [5]. This 
need led to the development of the Microstrip SQUID 
Amplifier (MSA) in which the input coil deposited on 
(but insulated from) the washer of a SQUID acts as a 
resonant microstrip [6]. Such amplifiers have achieved a 
noise temperature within a factor of two of the standard 
quantum limit [7-9]. More recently, new designs have 
appeared intended to extend the frequency of operation 
to frequencies as high as 10 GHz, aimed at the readout 
of superconducting qubits [10] and the detection of mi- 
cromechanical motion [11]. These include incorporation 
of a gradiometric SQUID at the end of a quarter wave 
resonator [12] and the direct injection of the microwave 
signal from a quarter wave resonator into one arm of the 
SQUID ring [13, 14]. 

Besides having desirable properties such as high gain, 
wide bandwidth and near quantum- limited operation, 
microwave SQUID amplifiers (MWSAs) — unlike conven- 
tional Josephson parametric amplifiers [15-17] — also of- 
fer an intrinsic separation of input and output chan- 
nels of the signal that makes them unique among am- 
plifiers based on Josephson tunnel junctions. This prop- 
erty makes them especially well suited as a preamplifier 
in the measurement chain for superconducting devices 
by eliminating the need for channel separation devices, 
such as circulators and isolators, between the sample un- 
der test and the first amplification stage. Although mi- 



crowave SQUID amplifiers have been successfully used 
experimentally, questions pertaining to their nonlinear 
dynamics and ultimate sensitivity as amplifiers have con- 
tinued to remain challenging problems. Previous the- 
ories include quantum Langevin simulations [4, 18] and 
treatment of the SQUID as an interacting quantum point 
contact [19]. The ultimate exploitation of the amplifier, 
however, requires a deeper understanding of its behavior 
at the Josephson frequency and its harmonics. Besides 
being valuable for practical considerations, such under- 
standing may help discern the cause of intrinsically non- 
reciprocal operation of the MWSA that has hitherto re- 
mained an open question. This concern is especially rel- 
evant to applications such as qubit readout where the 
amplifier backaction may prove to be the Achilles' heel. 
In this work, we develop an ab-initio theoretical frame- 
work to understand the high frequency dynamics of the 
SQUID in detail. In addition to giving us crucial in- 
sights into the amplifying mechanism of the MWSA and 
its nonreciprocal response between the input and output 
signal channels, this approach enables us to calculate the 
experimentally relevant quantities such as available gain, 
added noise and directionality at operating frequencies 
of interest. 

We perform our analysis in the paradigm of input- 
output theory and employ the method of harmonic bal- 
ance to study the driven dynamics of the device. The 
dc SQUID is biased in the voltage regime — in contrast 
to the usual Josephson parametric amplifiers operated in 
the zero voltage state with the phase excursions of the 
Josephson junction confined to a single cosine well — and 
has the dynamics of a particle sampling various wells of a 
two-dimensional tilted washboard [1]. The input-output 
analysis thus needs to be generalized to take into ac- 
count phase running evolution in this two-dimensional 
potential. Our approach involves a self-consistent deter- 
mination of the working point of the device established 
by static bias parameters (the static bias current Ib and 
external flux <£ ext shown in Fig. 1) followed by a study 
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FIG. 1: Circuit schematic of a conventional MWSA. The 
SQUID consists of two Josephson junctions, arranged in a 
superconducting loop, with inductance L. The loop is biased 
using a static current source 1b and an external flux ^ext- 
An input voltage V\ generates an oscillating current in an 
input coil inductively coupled to the SQUID thus inducing a 
small flux modulation 5& of the flux enclosed by the loop. For 
optimal flux bias [4> ex t = (2n + l)$o/4] that maximizes the 
flux-to- voltage transfer coefficient Vb = (dV/d& ex t)i B , this 
causes an output voltage Vi = V<f>5<& to develop across the 
ring. Thus the device behaves as a low impedance voltage 
amplifier. 



of the rf dynamics using a perturbative series expansion 
around this working point. In Sec. II we introduce our 
input-output model for the dc SQUID. Using this in Sec. 
Ill, we first derive the response at zero frequency and 
at the Josephson oscillation frequency uoj in a self con- 
sistent manner. Following this, we evaluate the pertur- 
bative response at finite frequency around zero and ujj 
as a scattering matrix in the basis of relevant modes of 
the circuit, which clearly elucidates the nonreciprocal dy- 
namics of the device. In Sees. IV and V we calculate the 
figures of merit such as power gain, directionality and 
noise temperature, and identify the fundamental limits 
on the performance of the device. Section VI contains 
our concluding remarks. 



II. ANALYTICAL MODEL 

A. SQUID circuit basics 

The SQUID circuit considered in our analysis is shown 
in Fig. 2(a). The dynamics of the system are modelled 
as a particle moving in a two-dimensional potential of the 
form [1] 
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Here, 1b is the bias current, Iq is the critical current 
of each junction, <^ ext = 27r^> ext /3>o represents the ex- 
ternally imposed flux in the loop, /3l = 2L1q/^q de- 
notes a dimensionless parametrization of the SQUID 



loop inductance, Ej = Io$o/2ir is the Josephson en- 
ergy and $0 = h/2e is the flux quantum. We have 
introduced the common, cp c = (ipL + ipn)/2, and dif- 
ferential, ip D = ((fL — ipr)/2, mode combinations of the 
phases of the two junctions that form the axes of the 
two-dimensional orthogonal coordinate system. 

To facilitate an input-output analysis of the circuit, we 
replace the resistive shunts across the junction with semi- 
infinite transmission lines [cf. Fig. 2(b)] of characteristic 
impedance Zq = R, following the Nyquist model of dis- 
sipation. Thus the shunts play the dual role of dissipa- 
tion and ports (or channels) used to address the device. 
This allows us to switch from a standing mode repre- 
sentation in terms of lumped element quantities such as 
voltages and currents to a propagating wave description 
in terms of signal waves travelling on the transmission 
lines. The amplitude of these waves is given by the well 
known input-output relation [20], 
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where V 1 and P denote the voltage across the shunt re- 
sistance and current flowing in the shunt resistance re- 
spectively. From Eq. (1), we obtain the common mode 
current, I c = (I L +I R )/2, and differential mode current, 
I D = (I L — l R )/2, flowing in the shunts by identifying 
(p c,D as the relevant position variables. The current in 
each mode can thus be interpreted as the "force" [21] 
that follows directly from Hamilton's equation of motion 
as 
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Here we have expressed the currents in equivalent fre- 
quency units, 
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and ujq = , (characteristic currents) 

(6) 



with cpo = ^>o/(27r). Including a capacitance across the 
junction gives an additional term, involving a second- 
order derivative of the common and differential mode 
fluxes, of the form —ujb^c^P C ' D with ft c = 2ttIbR 2 C /<&o, 
on the right-hand side of Eq. (4). This parametriza- 
tion of capacitance, motivated by calculational simplic- 
ity, leads to a different parametrization of the plasma 
frequency, uj p = (7 M)^) 1/2 [1 - (I B / Io) 2 ] 1/A ■ The more 
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FIG. 2: Equivalent input-output model of the SQUID (a) The bare SQUID loop without the input coupling circuit. As noted 
in Fig. 1, the circuit has two static biases - a common mode bias current Ib and a differential mode external flux $ ex t- There 
is also a capacitance C across each junction, not shown here for simplicity, (b) Equivalent SQUID circuit under Nyquist 
representation of shunt resistances and separate static current biases Ig and I§ for the left and right junction respectively. The 
common mode bias current Ib now corresponds to the even combination of two external bias currents — (Ib + while the 

external flux $ ex t corresponds to the differential combination of the two current sources L(l^ — Ib)/2- The oscillating signals 
are modelled as incoming and outgoing waves travelling on semi-infinite transmission lines, representing the shunt resistances 
across the two junctions, (c) Effective junction representation for evaluating the signal response of the device. Here, we have 
replaced the junctions biased with a static current with effective junctions pumped by the Josephson harmonics (represented 
by a "glowing" cross with a pumping wave) generated by phase running evolution in the voltage state of the junction. 



conventional parametrization with a fixed value of ca- 
pacitance for all bias values can be implemented in a 
more comprehensive calculation aided by numerical tech- 
niques. 

Equation (4) represents a subtle current-phase rela- 
tionship for the two-junction system, analogous to the 
first Josephson relation. We note that Eqs. (4) can al- 
ternatively be derived using a first-principles KirchofT law 
analysis of the circuit in Fig. 2(a). Similar to the cur- 
rents, we can define the common and differential mode 
voltages as 

£ c = — ; S^eL. (7) 

Further, by the second Josephson relation, we have 

(ou c ) = V dc /ipo = ujj, (8) 

where Vdc is the static voltage developed across the 
SQUID biased in the running state. 

We note that the usual mode of operation of a dc 
SQUID involves an input flux inductively coupled using 
an input transformer of which the loop inductance forms 



the secondary coil (Fig. 1). The input transformer, how- 
ever, is an experimental artifact required to ensure the 
impedance matching with the input impedance of the 
SQUID at a desired frequency. It is not crucial from the 
point of view of device characteristics, however, as it is 
the SQUID which provides amplification and all the rel- 
evant nonlinear dynamics of the device. In the ensuing 
analysis we do not employ a separate input port, but 
rather consider a direct input coupling through the dif- 
ferential mode of the ring which couples to the flux in 
an analogous manner [Fig. 2(b)]. Such a scheme may 
also prove beneficial for a practical device to overcome 
the problem of low coupling at high signal frequencies, 
as recently shown experimentally using a SLUG (Super- 
conducting Low-inductance Undulatory Galvanometer) 
microwave amplifier [13, 14]. 

B. Harmonic balance treatment 

Using the input-output relation of Eq. (2) with Eqs. 
(4) and (7), we obtain the equations 

u C ' D (t) = Q C ' D (t) + 2uj inC > D (t) (9) 
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FIG. 3: (Color online) Spectral density landscape of com- 
mon and differential modes of the SQUID. The tall solid ar- 
rows show the Josephson harmonics generated internally in 
the running state of the device. The small input signal fre- 
quency com and different sidebands generated by mixing with 
Josephson harmonics are shown with dashed arrows. 



for common and differential mode circuit quantities. We 
employ the technique of harmonic balance and solve Eq. 
(9) in the frequency domain, at all frequencies of interest 
(see Fig. 3). This is achieved by assuming two parts to 
the solution for each variable of interest (ip c and (p D ), 



if D = <fo + 5tp D (t), 



(10) 
(11) 



where uojt and 0o represent the average static values of 
the common and differential mode phases [cf. Eq. (8)]. 
The time varying components are of the form 
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where H(t) refers to the components at the Josephson 
frequency ujj and its harmonics. The term £(t) includes 
the components oscillating at the signal frequency uj m 
and its resultant sidebands uo n = nuoj + uj m generated by 
wave mixing via the nonlinearity of the SQUID: 
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We note that the number of Josephson harmonics in- 
cluded in the analysis [i.e. K in Eq. (13)] is determined 
by the order of expansion of the junction nonlinearity 
in Sip. This in turn is determined by the bias voltage 
of the device set by the bias current Ib- As Ib is re- 
duced towards the critical current of the junction Jo, 
higher Josephson harmonics become more significant as 
the characteristics of the device become increasingly non- 
linear. We can, therefore, calculate the response pertur- 
batively by expanding each of the coefficients p and s in 
Eqs. (13)- (14) as a truncated power series in the reduced 
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The degree of the resultant polynomial evaluation of p, s 
coefficients is set by the desired order of expansion in Sip. 
As e < 0.5 (or equivalently I B > 2I ) for the SQUID 
to operate in the running state at any value of flux bias 
[1], which is the regime of interest for the SQUID to be 
operated as a voltage amplifier, it provides a convenient 
small parameter of choice ensuring rapid convergence of 
the perturbation series method. Furthermore, this pa- 
rameter serves as the effective strength of the different 
Josephson harmonics which play a role analogous to the 
strong "pump" tone of conventional parametric ampli- 
fiers. 



III. CALCULATION OF SQUID DYNAMICS 

A. Steady state response: I — V characteristics 

We first determine the working point of the SQUID by 
solving for the steady state characteristics. As the zero 
frequency response of the system is intrinsically related to 
the response at the Josephson frequency through Eq. (8), 
we calculate it self-consistently along with the strength 
of the various Josephson harmonics in the steady state 
by considering only the static source terms with no os- 
cillating input drive at the Josephson frequency and its 
harmonics. This yields a set of boundary conditions of 
the form 



)[kuj] - ud[kuoj} = k G [0, K]. 
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We solve this set of simultaneous equations to calculate 
the strength of the various Josephson harmonics gener- 
ated internally from the static bias due to the junction 
nonlinearity along with the zero frequency characteris- 
tics. Figure 4 shows a plot of the static transfer function 
= d{uo c ) /dip ex t obtained using the perturbative se- 
ries method to determine the coefficients p^ D [Eq. (13)] 
described in the last section. The agreement between 
the exact numerical calculation and the perturbative an- 
alytical calculation improves on increasing the order of 
the perturbation series expansion by including mixing 
processes mediated by higher Josephson harmonics. Fur- 
ther, from the steady state calculation for the differential 
mode, we obtain a relation for the phase angle between 
the two junctions in the ring as 
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where the coefficients a& are of order unity. Thus, we see 
that the average values of both the explicit static bias pa- 
rameters namely e (common) and cp ext (differential) par- 
ticipate in establishing each of the implicit static biases 



4 




0.1 0.2 0.3 0.4 0.5 0.6 

8 



0.7 




FIG. 4: (Color online) Static transfer function of the SQUID 
calculated as a function of two bias parameters e = ujq/ujb 
at ifext = tv/2 for (a) strongly overdamped (Qc = 0) and 
(b) intermediately damped junctions (Qc = 1)- The (black) 
triangles represent the transfer function calculated from the 
exact numerical integration of the SQUID equations. The 
(green) circles correspond to the K = 1 evaluation includ- 
ing only the Josephson frequency [Eq. (13)]. This first order 
evaluation does not show any voltage modulation with flux 
as there is no coupling between the common and differential 
modes at this order. The (blue) squares and (red) diamonds 
correspond to an evaluation including the second (K = 2) and 
third Josephson (K — 3) harmonic respectively. The corre- 
sponding curves represent interpolating polynomials. In both 
plots, the agreement of the perturbative series with the exact 
numerical solution improves on including higher order cor- 
rections corresponding to contributions of higher Josephson 
harmonics. 



- Vdc (or equivalently ujj) for the common mode and </>o 
for the differential mode. The contributions arising from 
the bias current, as shown in Eq. (17), lead to a rolling of 
the static phase difference around the SQUID loop that 
manifests itself as the change in curvature of the transfer 
function curves shown in Fig. (4). Furthermore, we note 
that, as indicated by the steady state calculation, the 
flux dynamics of v$ evaluated using the truncated har- 
monic series calculation are 'slower', that is, they shift 
to higher values of bias with respect to the exact numer- 
ical results. Nonetheless, the predicted magnitudes are 
comparable and hence the theory is capable of making 
semi-quantitative predictions in an analytically tractable 



manner. The major merit of this approach over conven- 
tional methods lies in the natural extension offered for 
the study of higher frequency dynamics as discussed in 
the following sections. 

B. RF response: Scattering Matrix 

Once we have determined the static working point for 
the SQUID, we can solve for its rf dynamics in the small 
signal regime. The aim is to calculate signal amplitudes 
by including the Yi C ' D (t) term in our analysis and con- 
sidering all the mixing processes mediated by the pumps 
U CjD (t) evaluated in the last section, permissible by the 
harmonic balance of Eqs. (4), (7) and (9). This is equiv- 
alent to the representation shown in Fig. 2(c), where we 
model the mixing of the input signal by the SQUID as a 
parametric interaction with different Josephson harmon- 
ics playing the role of an effective "colored" pump. In 
the limit of a small amplitude input signal, which is the 
relevant limit for most practical situations, we can then 
introduce a linear response description of the dynamics 
as an admittance matrix seen from the ports. This can 
be obtained from the current-phase and voltage-phase 
relationship [Eqs. (4), (7)] as 

~u = (18) 

and 

= Mjt, (19) 

yielding 

Y = MM -1 . (20) 

The vectors in the equations above are defined in the 
basis of all signal and sideband frequencies of interest, 
(Z c [nujj + cj m ], Z D [nujj + u m ]), ne [-N, leading 
to a 4(27V + 1) x4(27V + l) admittance matrix. We further 
note that the matrix is block diagonal since harmonic 
balance leads to two disjoint manifolds, each of which 
forms a closed subspace of dimension 2(2iV + 1). 

From the admittance matrix of Eq. (20), we can evalu- 
ate the scattering matrix of the SQUID using the identity 

S = (U + Y) _1 (U - Y). (21) 

Figure 5 shows the calculation for different orders in junc- 
tion nonlinearity and the relevant forward (\s CD \ 2 ) and 
backward scattering gain (l^^^ | 2 ). It immediately shows 
the emergence of the nonreciprocal gain of the device 
that, unlike conventional par amps, enables a two-port 
operation. As the nonlinearity of the device character- 
istics is increased by reducing 1b towards Io (thus in- 
creasing the expansion parameter e), we need to include 
the higher Josephson harmonics in the calculation which 
become significant due to rapid running evolution of the 
phase of the junctions in the two-dimensional tilted wash- 
board. This leads to a situation analogous to pumping of 
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FIG. 5: Josephson harmonics and small signal scattering gain of the SQUID calculated using harmonic balance with expansion 
of the sin<£ nonlinearity to (a) first, (b) third and (c) fifth order respectively. Each panel shows the relevant modes of the 
frequency spectrum included in the calculation at that order [Eqs. (13)-(14)]. The dispersive mixing between various temporal 
modes of the system is denoted using grey arcs with the relevant Josephson harmonic acting as the pump indicated next to 
them. The relative strength of the different mixing processes is indicated by the respective widths of the arcs, with the strongest 
being denoted by the thickest arcs. Also shown are plots of II(£), the effective pumps in common (blue) and differential (red) 



modes at each order of the calculation. The box panels show the respective forward (\s 



and backward (|s 



scattering 



gains as a function of reduced input frequency cj m /cJo and bias parameter e — uoq/uob- The surface plot in (a), calculated 
using only the Josephson frequency, shows no asymmetry between the forward and backward gains (blue and red surface 
plots respectively) . The asymmetry develops on inclusion of higher harmonics that implement a multitone pump which is not 
symmetric about t = 0, as seen from the plots of IJ(t) in panels (b) and (c). As we increase the order of calculation and include 
higher harmonics, the asymmetry increases and finally peaks at an optimal value of bias parameter e = 0.455. 



the SQUID by an effective multitone pump of the form 
II(t) = Y^k=i P k cos (ujt + fa) in both C and D modes 
[see Il(t) panels in Fig. 5]. The dynamics of such a 
system include multi-path interference involving differ- 
ent Josephson harmonics. This effect, analogous to sym- 
metry breaking in ratchet physics [22], implements an 



asymmetric frequency conversion scheme guided by rela- 
tive phases <pk of different Josephson harmonics driving 
the junctions [23]. The signal in the differential mode is 
preferentially upconverted, coupled through higher order 
mixing processes into the common mode and then pref- 
erentially downconverted into the common mode, yield- 
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ing a net forward gain from the differential mode to the 
common mode. The reverse gain process from C to D is 
disfavored by the same reasoning, leading to the nonre- 
ciprocal operation of the SQUID amplifier. 



quantities of Eq. (22). Using the Y- matrix, derived in 
Eq. (20), we can write the impedance matrix Z of the dc 
SQUID as 



(U + Y) _1 



(24) 



IV. POWER GAIN OF THE SQUID 

The dc SQUID operated as a two-port voltage am- 
plifier resembles the configuration of a semiconductor, 
operational amplifier (op-amp) as opposed to that of a 
conventional parametric amplifier, which is a matched 
device. (That is, the input and output impedances are 
identical to the impedances of the transmission lines or 
coaxial cables). In this sense the MWSA is the magnetic 
dual of the rf SET (single electron transistor) [24]. The 
dc SQUID amplifies an input current (directly coupled as 
in this analysis or coupled as a flux via an input trans- 
former), and has a much lower impedance than the elec- 
tromagnetic environment in which it is embedded. Con- 
versely, the rf SET amplifies an input voltage, and has 
a much higher impedance than the electromagnetic envi- 
ronment in which it is embedded. The true power gain of 
either device, as seen from the ports, thus involves a de- 
embedding of the device characteristics. In the case of the 
SQUID, this requires a translation from the matched (or 
scattering) description based on the input-output theory 
considered in this paper to the op-amp or hybrid repre- 
sentation that is well suited for describing an unmatched 
amplifier such as the microwave MWSA. 

The hybrid matrix describing a two-port amplifier is 
of the form [25] 



with 
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where (Vi, I\) and (V2, h) denote the voltage and current 
associated with the input and output ports respectively. 
The power gain for such an amplifier is given by 
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where Ay is the voltage gain of the amplifier, Y in is the in- 
put admittance and Z out is the output impedance. Equa- 
tion (23) represents the gain of an effective "matched" 
device accounting for the impedance mismatch at the in- 
put and output ports. 

In principle, although the calculation of quantities in 
Eq. (22) can be performed using the scattering matrix 
evaluated in Eq. (21) [25], nonetheless it is advantageous 
to transform to a description that is more natural in de- 
scribing the relationship between standing mode current 
and voltage variables. We find that an impedance matrix 
(Z) representation is well suited for such a purpose due to 
its rather straightforward mapping to the standing mode 




(25) 



Here, as before, uj , Q are vectors defined in the space 
of all signal and sideband frequencies of interest. Also U 
is an identity matrix of appropriate dimensions and cor- 
responds to the admittance contribution of the resistive 
shunts across the junctions. 

The next step is to make the translation from the 
impedance matrix derived in the common and differential 
mode basis to the two-port description of Eq. (22). This 
requires an identification of the correct "input" and "out- 
put" voltages and currents for the circuit in Fig. 2(a). As 
the SQUID readout involves measurement of the voltage 
developed across it, the relevant output quantities are re- 
lated to the common mode quantities as V2 = V c and 
I2 = 2I C . The translation to the input variables of the 
hybrid representation is more subtle. For this purpose 
we first note that, in conventional SQUID operation, the 
input flux coupled into the ring modulates the circulating 
current J which is, thus, the relevant input current of the 
device. The equivalent input voltage that causes the flux 
modulation of the circulating current can be represented 
by a voltage source Vj in series with the inductance of 
the loop. Figure 6 summarizes the different possible two- 
port representations of the SQUID used in this paper. 

On interpreting the loop variables (Vj, J) described 
above in terms of the differential mode voltage V D and 
current I D (see appendix A for details), we obtain the 
following equivalence between the coefficients of the hy- 
brid matrix in Eq. (22) and the Z- matrix of Eq. (25): 
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(27) 
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Using the above translation in Eq. (23), we find an ex- 
pression for the power gain purely in terms of Z-matrix 
coefficients: 



Gp [id m ] 



^CD 



^m]| 2 



Re\z cc \(jj m me\z DD \u> n 



(30) 



Figure 7 shows the power gain of the device as a function 
of bias and input frequency, calculated using Eq. (30). It 
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H W 2 §\i 2 



X 



Vj/2 Vj/2 




FIG. 6: Different representations of a two-port network and 
analog configurations for the dc SQUID, (a) Y-matrix rep- 
resentation defined for closed boundary conditions Yij = 
dIi/dVj\v k ^j=o for the junctions and inductance, omitting 
the shunts [Eq. (20)]. (b) Z-matrix representation defined 
for open boundary conditions Zij = dVi/dIj\i k , j= o including 
the shunts [Eq. (24)]. (c) (Hybrid) H-matrix or op-amp repre- 
sentation defined with mixed boundary conditions [Eq. (22)]. 
In effective matrices for the SQUID, the common mode (C) 
and differential mode (D) excitations of the ring play the role 
of ports 1 and 2, if the SQUID is addressed using hybrids. In 
each panel, the quantities shown with solid arrows represent 
the stimulus while those shown with dashed arrows represent 
the corresponding response of the network. 



shows that power gain of the MWSA increases quadrati- 
cally with decreasing input signal frequency, a result cor- 
roborated by a simple quasistatic treatment presented 
in Appendix B. The reverse power gain of the device is 
calculated in a similar manner as 
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FIG. 7: Power gain of the MWSA calculated with K = 3, 
N = 2, taking into account the modification of the input and 
output impedances of the device by matched loads. The pa- 
rameters are <I> = <E>o/4, 0l = 1 and Qc = 1. (a) Power 
gain versus bias parameter s = ujq/ujb calculated for a fixed 
input frequency cj m = 0.01 ujq. The solid curve is an in- 
terpolating polynomial of degree two. (b) Power gain versus 
input frequency uj m /oJo calculated with bias parameter fixed 
at £ = 0.455, the optimum value for attaining minimum noise 
temperature [see Fig. 9(a)] at low frequencies (u <C wo). The 
fit is of the form Gp = [0.006/ (u m /wo) 2 ] + 2 measured in 
linear units. 



i?Re[Z h 



A? 



J|Re[Z, 



out J 
DC\ 



Re[F in ]Re[Z out ] 

, 112 



[Eq. (22)] 



Re[z 



cc 



UJ ri 



}Re[z 



DD 



(31) 



The directionality (Gp — G r p v ) — which is a measure of 
the asymmetry between forward and reverse power gains 
— follows directly from the asymmetric scattering gain 
discussed in the previous section. Our calculation shows 
that it is a strong function of the bias e (Fig. 8); further- 
more the optimal bias for maximum power gain is not the 
same as that for maximal directionality. We note that the 
results presented here have been obtained with a trun- 
cated harmonic series excluding all Josephson harmonics 
above 3ujj. In the real device, the achievable isolation 
between forward (differential-to-common) and backward 
(common-to-differential) gain channels may be quantita- 
tively different due to the presence of the neglected higher 
order interferences. 



V. NOISE TEMPERATURE 

In this section, we evaluate the noise added by the dc 
SQUID operated as a voltage amplifier. The noise added 
by a system can be quantified by its noise temperature, 
T^, defined as 



k B 



(32) 



where A is the Caves added noise number [7]. This noise 
temperature corresponds to the effective input temper- 
ature of the amplifier obtained by referring the added 
noise measured at the output to the input, and is quan- 
tified in terms of energy quanta per photon at the sig- 
nal frequency. For a phase preserving amplifier, such 
as the MWSA, the minimum possible noise temperature 
corresponds to half a photon of added noise, that is, 
^-min = 0.5 in the large gain limit (in general, A m i n = 
1/2- 1/2G). 



8 



0.00 



35 dB 
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FIG. 8: Directionality (in dB) of the MWSA as a function 
of bias parameter e — ujq/ujb and reduced input frequency 
oJm/oJo- The red dots represent high directionality and blue 
dots represent low directionality. The parameters used were 
the same as those in Fig. 7. 



Using the hybrid representation developed in the previ- 
ous section and Appendix A, we write the noise inequality 
for the MWSA as 



k B T N > 



ycycbjj 



Re[S v cj]i-Im[S v cj] 



, (33) 



where Syv represents the spectral density of the voltage 
fluctuations at the output, Sjj represents the spectral 
density of the circulating current fluctuations and Syj 
is the cross-correlation between the voltage and current 
fluctuations [25]. 

As in the case of power gain, we can evaluate these 
quantities from the Z matrix of the SQUID derived in 
Sec. IV. This exercise is enabled by the fact that the 
input-output theory treats the deterministic signal in- 
put and noise of the system on an equal footing. Thus, 
the linear response description developed to calculate the 
signal gain provides a straightforward way of generaliz- 
ing the theory to understand the noise properties of the 
system, simply by replacing the input current signal with 
a noise signal described by a spectral density of the form 



S n [u>] = 2hujRe[Y[uj}} coth (^^) 



(34) 



As before, we use the Z matrix to calculate the spectral 
densities in Eq. (33). Using this we write the voltage 
noise spectral density in the common mode as 

N 

Sycyc = ^ \Zon\ 2 Sl c I c [ nUJ J+Um] 
n=-N 

N 

+ E \zon D \ 2 Si»i°[nuj + ou m }. (35) 

n— — N 

Here, the first sum accounts for the contribution to the 
noise arising from the common mode signal (n = 0) and 



sidebands about the Josephson harmonics (n = ±1,±2) 
included in the calculation; the second sum accounts for 
the noise generated in the common mode output signal by 
the differential mode signal and sidebands, arising from 
coupling between C and D modes. Similarly, we can 
calculate Sjj as 



S. 



j j 



JyDyD , 



(36) 



by making the identification J = 2V D / \iujL). Here, as 
before, we calculate SyD V D from the Z matrix: 



N 



S V d v d = ^ IZq^^SjCjcIuUJj + 0^ 



-N 



N 



+ E \Zon D \ 2 Sl*in[™J+UJ m ]. (37) 



-N 



Finally, for Syj we have 
/ n 



*>v c J 



J2 Z0nZ0n*SlClc[nUJj+UJ„ 



N 



-N 



+ E zZ D Zon D *S I »Mnwj + u m ] .(38) 

n= — N ) 

Figure 9 shows plots of the Caves noise number of the 
device, calculated in both the thermal regime [fc#T ^> 
ftijj m , where all terms in Eqs. (35)- (38) contribute equal 
noise powers] and the quantum regime [fc#T <C ^cj m , 
where each term in Eqs. (35)- (38) contributes a noise 
power proportional to its frequency in accordance with 
Eq. (34)]. 

There are a number of points to be highlighted. Our 
calculation shows that, at the optimal current bias of 
e = 0.455, the MWSA attains the quantum limit of 
added noise corresponding to half-photon at signal fre- 
quency. Moreover, the optimum bias point for minimum 
noise corresponds to the bias for maximum scattering 
gain [Fig. 5(c)] rather than for the maximum power gain 
[Fig. 7(a)]. This result, previously found both theo- 
retically and experimentally [26], follows from the fact 
that the added noise is a property of the bare SQUID 
without any matching to input and output loads. In 
the case of conventional parametric amplifiers, the min- 
imum noise indeed occurs at the maximum scattering 
gain. Furthermore, the partial cross-correlation between 
the output voltage noise across the SQUID and the su- 
percurrent noise circulating in the loop is crucial to min- 
imizing the noise in both thermal and quantum regimes. 
We also note that, for sufficiently low signal frequencies, 
the calculated added noise number is found to saturate 
at a value slightly below the quantum limit of one half- 
photon at the signal frequency. This result, we suspect, 
is due to the fact that at the bias for minimum noise, 
the reverse gain is substantial and hence the isolation is 
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FIG. 9: (Color online) Caves added noise number for the 
MWSA, calculated using the harmonic balance analysis with 
K — 3, N — 2 as a function of (a) e — ujq/ujb for 
ujm = 0.01 ojo and (b) reduced input frequency uj m /oJo with 
e — 0.455. In both plots, (round) black markers show the 
noise number A = T/v/T obtained in the thermal regime 
ksT ^> hum, while (square) green markers show the noise 
number A = JzbTn /hujm calculated in the quantum regime 
ksT <C hujm- The solid curves in (a) represent interpolat- 
ing polynomials. The quantum calculation gives a minimum 
value for A = JzbTn /hujm ~ 0.5, attained at e = 0.455, corre- 
sponding to one-half photon of added noise [horizontal black 
dashed line in (a)]. The optimal value of bias current for min- 
imum added noise does not coincide with that for achieving 
the maximum power gain [Fig. 7(a)] or directionality (Fig. 



not perfect (Fig. 8). The quantum limit of one half- 
photon is a limiting value calculated for ideal detectors 
with zero reverse gain and high forward gain [25], a con- 
dition which is not satisfied at the optimal noise bias in 
our calculation. Finally, our calculation shows that the 
minimum noise number is achieved only when the signal 
frequency is much lower than the characteristic Joseph- 
son frequency ujo = 2ttIoR/ <£>o, and increases significantly 
with increasing signal frequency. 



VI. CONCLUDING REMARKS 

In summary, we have developed a new method based 
on input-output theory to provide a first-principles anal- 



ysis of the microwave SQUID amplifier (MWSA). In this 
paradigm we treat the SQUID biased in its running state 
as a parametric amplifier pumped by a combination of 
Josephson harmonics generated internally by the motion 
of the phase of the junctions. This approach leads to a 
fully self-consistent description of both the static and rf 
dynamics of the device. The scattering matrix calcula- 
tion shows that the nonreciprocal gain of the amplifier 
arises from mixing processes involving higher Josephson 
harmonics which implement an asymmetric frequency 
conversion scheme involving upconversion to and subse- 
quent downconversion from the Josephson frequency. We 
find that the power gain of the matched SQUID ampli- 
fier decreases quadratically with signal frequency u m ; by 
comparison, the gain in the usual SQUID operation with 
a matched input coil scales as l/uo m [26]. However, a re- 
cently reported dc SQUID amplifier [14] using the direct 
coupling method considered in this paper demonstrated 
that the power gain scaled as at a frequency of a 

few GHz and a bandwidth of several hundred MHz. 

Our analysis shows that the MWSA achieves quantum- 
limited noise performance for optimal flux and current 
biases, and for signal frequencies significantly lower than 
the characteristic Josephson frequency uoq = 2ttIqR/^q. 
The added noise increases significantly with increasing 
frequency. This problem can be alleviated by using junc- 
tions with higher values of critical currents. With the 
present technology for niobium junctions, critical current 
densities of tens of microamperes per square micron are 
readily achievable. This translates into characteristic fre- 
quencies of about 100 GHz, which should be sufficient to 
achieve lower noise at GHz frequencies provided hot elec- 
tron effects due to dissipation in the shunts are mitigated 
[27]. Furthermore, our analysis shows that simultaneous 
optimization of gain, directionality and noise is a delicate 
operation since the optimal biases for these three prop- 
erties do not coincide. Based on our calculation, at the 
working point for minimum added noise, A « 0.5, power 
gains of 15 — 18 dB and directionality of around 5 — 8 dB 
are obtained. However, higher power gains of 20 — 30 dB 
and directionality of 10 — 12 dB can be realized by per- 
mitting a higher noise number A « 5 — 10. Though the 
predicted directionality is still modest, it suffices to re- 
duce the number of nonreciprocal elements (circulators, 
isolators) in the measurement chain typically employed 
for the readout of superconducting qubits. Moreover, the 
noise penalty incurred with MWSAs compares very well 
to standard cryogenic amplifiers such as HEMTs whose 
typical noise numbers lie in the range 40 — 50 for mi- 
crowave frequencies. 

Although the results presented in this paper are semi- 
quantitative we believe that extension of the analysis to 
higher orders, in conjunction with numerical optimiza- 
tion techniques, can be a useful tool to analyze SQUID- 
based devices due to rapid convergence offered by the 
harmonic series method. This approach would allow one 
to evaluate the appropriate parameters, depending on 
the intended application, that yield the best compromise 
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between gain and noise properties. 
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Figure 10 shows the two representations [Figs. 6(b) 
and (c)], one in terms of differential mode quantities 
(V D , I D ) suitable for a scattering or matched representa- 
tion (since the input and output impedances are just the 
transmission line impedance) and the other in terms of a 
circulating current J and a loop voltage Vj , which are the 
relevant input quantities for the device in an unmatched 
hybrid description. In Fig. 10(a) Kirchoff's current law 
gives 

J = I d -Il, (Al) 
while in Fig. 10(b), from Kirchoff's voltage law, we have 



V 



-J, 



(A2) 



Appendix A: Loop Variables for the dc SQUID 

In this appendix we establish the correspondence be- 
tween the differential mode variables that serve as the 
input in the analysis of Sees. II and III and the input 
variables required for the hybrid representation discussed 
in Sec. IV. The output variables in the two representa- 
tions have a simple relationship as explained in Sec. IV. 



v=o 




FIG. 10: (Color online) Equivalence between the SQUID dif- 
ferential mode and the op-amp input variables. Here, for 
simplicity, we have used the symmetric version of the SQUID 
to divide the ring along the equipotential. The circuit in (a) 
models the device as an impedance response function to an 
imposed current source I D , as a result of which a voltage drop 
+V D (-V D ) develops across the left (right) junction, (b) Hy- 
brid representation for the SQUID, which models the input 
response by introducing a differential voltage source in the 
SQUID loop and recording the current that flows through the 
junction. The junction at this point is replaced with an effec- 
tive junction pumped using the various Josephson harmonics 
generated by the static bias current [also see Fig. 2(c)]. 



with J = 1^ . 

To establish the equivalence of the two representations 
from the point of view of the junction, we require the 
voltage across the junction V D and current through the 
junction J to be conserved (see Fig. 10). Thus, using 
Eq. (Al) in Eq. (A2), we obtain 



Vj = -iuLI D . 



(A3) 



Similarly it is easily seen that the circulating current J 
is given as 

2V D 



J 



iujL 



(A4) 



Appendix B: Static analog circuit for the SQUID 

The SQUID can be thought of as a current amplifier 
with a current transferred from a low-impedance input 
port to a high-impedance output port. This description 
is analogous to the FET dual model with the gain given 
by a transimpedance instead of a transconductance. The 
equivalent 'current gain' of such a device (Fig. 11) for 
frequencies sufficiently close to zero [cj m <C p in R/L = 
ojo/(ttPl) to be precise] can be modelled as 



-'out 



~Rd' 



This leads to a power gain 



G d P c = 



Rd 



Re[Z h 



(Bl) 



(B2) 



For frequencies of interest, Re[Zi n ] w uj^L 2 / (p- in R) . Us- 
ing this result in Eq. (B2), we obtain the power gain 



G 



dc 



Pin / Vj> 
Pout V^ra 



which can be rewritten as, 



cup 

UJ 1 



(B3) 



(B4) 
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Differential Common 



FIG. 11: Equivalent low frequency circuit for a SQUID for cal- 
culation of unilateral power gain. The input circuit is mod- 
elled as an effective impedance viewed by a low frequency 
differential mode current. The output circuit impedance 
comprises a bias-dependent resistor, denoting the dynamic 
impedance of the junction, that converts the output voltage 
to a corresponding output current. The net "transimpedance" 
is given by the static flux-to- voltage transfer function of the 
device. The symbols pin, out denote bias-dependent constants 
of order unity. 

where p g is a bias-dependent and frequency-independent 
constant of order unity. Here, we have used the relation 
K Pt = R/L = ujq/tt [2] for f3 L = 1. Equation (B4) shows 
that the gain drops quadratically with increasing signal 
frequency, and that no power gain is obtained for signal 
frequencies close to the plasma frequency of each junction 
in the SQUID. This frequency dependence of the power 
gain is borne out by the full rf analysis shown in Fig. 
7(b). Figure 12 shows a comparison of gain calculated 
using quasistatic response functions as shown in Eq. (B3) 



and a rf calculation at low frequencies, involving the third 
Josephson harmonic [same as that shown in Fig 7(a)]. 
The agreement is better for lower values of e where high 
frequency components of the device are less significant. 
The impedance matrix calculation generates extra terms 
due to self-summation caused by the inversion operation 
[cf. Eq. (24)] which leads to higher order corrections 
absent from the purely quasistatic calculation. 
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FIG. 12: Comparison of the power gain as a function of bias 
e — ujq/ujb calculated using a K — 3, N — 2 calculation 
(solid black line) and a purely quasistatic calculation (dashed 
red line) of Eq. (B3). For the quasistatic gain calculation, 
Vb and pi n ,out were obtained from the I — V characteristics 
evaluated in Sec. Ill A. 
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